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Abstract. The surface temperature of the Greenland ice sheet 
is among the most important climate variables for assess- 
ing how climate change may impact human societies due 
to its association with sea level rise. However, the causes 
of multidecadal-to-centennial temperature changes in Green- 
land temperatures are not well understood, largely owing to 
short observational records. To examine these, we calculated 
the Greenland temperature anomalies (GTA[g-nhj) over the 
past 800 yr by subtracting the standardized northern hemi- 
spheric (NH) temperature from the standardized Greenland 
temperature. This decomposes the Greenland temperature 
variation into background climate (NH); polar amplification; 
and regional variability (GTA[g-a:h]). The central Green- 
land polar amplification factor as expressed by the variance 
ratio Greenland/NH is 2.6 over the past 161 yr, and 3. 3-4. 2 
over the past 800 yr. The GTAjg-nh] explains 31-35 % of the 
variation of Greenland temperature in the multidecadal-to- 
centennial time scale over the past 800 yr. We found that the 
GTA[g-nh] has been influenced by solar-induced changes in 
atmospheric circulation patterns such as those produced by 
the North Atlantic Oscillation/ Arctic Oscillation (NAO/ AO). 
Climate modeling and proxy temperature records indicate 
that the anomaly is also likely linked to solar-paced changes 
in the Atlantic meridional overturning circulation (AMOC) 
and associated changes in northward oceanic heat transport. 


1 Introduction 

The surface air temperature over Greenland is known to be 
affected by drivers such as large volcanic eruptions, changes 
in solar output, or anthropogenic forcings and by variability 
patterns such as the North Atlantic Oscillation (NAO) and/or 
the Atlantic multidecadal oscillation (AMO) (Hurrell, 1995; 
Hanna and Cappelen, 2003; Vinther et ah, 2003; Box et ah, 
2009). These regional variations have caused Greenland tem- 
perature to depart from the baseline trend defined by the av- 
erage northern hemispheric (NH) temperatures (Box et ah, 
2009). For more accurate projections of the future mass bal- 
ance of the Greenland ice sheet and of sea level changes, it is 
critical to understand the mechanisms underlying Greenland 
surface temperature changes. 

Greenland surface temperature, especially the temperature 
in western Greenland, fluctuates in a manner inverse to the 
observed NAO index (Hanna and Cappelen, 2003; Vinther et 
ah, 2003; Box et ah, 2009). The NAO, an important mode of 
climate variability in the North Atlantic that influences North 
Africa, Europe and North America, particularly in the win- 
ter, is often represented by the difference in sea level pres- 
sure between the Azores high and the Icelandic low (Hur- 
rell, 1995). Recent studies have shown that NAO variabil- 
ity is closely linked to changes in the atmospheric block- 
ing regime, which constitutes quasi-stationary anti-cyclonic 
systems interrupting the normal zonal flow (Barriopedro et 
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al., 2008; Woollings et al., 2010; Hakkinen et al., 2011; 
Rimbu and Lohmann, 2011). When the NAO index is posi- 
tive, west Greenland temperature in particular decreases due 
to stronger northerly winds, whereas northwestern Europe 
experiences warmer temperatures in response to stronger 
westerlies, and when the NAO index is negative, the oppo- 
site effects are observed (Wanner et al., 2001). An influence 
of the NAO on the average NH temperatures has also been 
suggested (Hurrell, 1996). Climate models (Shindell et al., 
2001) and observations (Kodera, 2002; Ineson et al., 2011) 
indicate a link between NAO-like variability and changes in 
solar output (Gray et al., 2010). 

During periods of stronger solar activity, ozone in the up- 
per stratosphere and lower mesosphere absorbs the increased 
ultraviolet irradiance, which increases the ambient tempera- 
ture. Higher ultraviolet irradiance also increases ozone pro- 
duction through the photodissociation of molecular oxygen 
in the stratosphere (Haigh, 1999), and the enhanced ozone in 
turn absorbs more radiation and further increases tempera- 
ture. Temperature changes maximize in the tropical latitudes, 
so that the temperature gradients between low and high lati- 
tudes in the stratosphere are amplified, creating stronger tro- 
pospheric westerly winds with associated pressure and tem- 
perature fields extending down to the surface where they 
resemble positive NAO-like atmospheric patterns (Shindell 
et al., 2001; Kodera and Kuroda, 2002). During periods of 
lower solar activity such as the Maunder Minimum, climate 
models and proxy records both suggest negative NAO-like 
patterns and warm southwest (especially southern edge of) 
Greenland temperature anomalies (Shindell et al., 2001). The 
proxy-derived spatial pattern of temperature change from the 
Medieval Warm Period to the Little Ice Age was shown to 
possess some features resembling the NAO/ AO (Arctic Os- 
cillation) pattern (Mann et al., 2009). 

In this study, we investigated the variations in the Green- 
land temperature anomaly using a recently developed Green- 
land temperature record (Kobashi et al., 2010, 2011) derived 
from trapped air in the GISP2 Greenland ice core from the 
Greenland Summit station and using two proxy average NH 
temperature reconstructions (Moberg et al., 2005; Mann et 
al., 2009). Proposed mechanisms for the temperature varia- 
tions were then investigated with a climate model (GISS-ER) 
and grid proxy reconstruction of global temperatures (Mann 
et al., 2009). We note that the inferences investigated in this 
paper were further tested and developed in a companion pa- 
per (Kobashi et al., 2012) dealing with the past 4000 yr of 
Greenland temperature record (Kobashi et al., 2011). 

2 Correlation coefficients and significance 

In this study, the Pearson product-moment correlation co- 
efficient was calculated to evaluate the strength of a corre- 
lation. To take into account the autocorrelation of time se- 
ries, the number of effective degrees of freedom (df e ) for p 


statistics was estimated from the effective sample size ( N t ), 
which depends on the effective decorrelation time (T c ) after 
Ito and Minobe (2010). The two relevant time series were 
represented by X and Y , and the sample size is N. Here, 
the time interval (At ) was 1 yr. Therefore, N e = N At / T e and 
df e = N t —2. For the correlation coefficients, T t was esti- 
N 

mated as T e = At J2 (Pxx(r) pyy(t) + Pxy(?) Pyx(?)), 
t=—N 

where r is a lag or a lead in years, and p is the correlation 
coefficient of the two time series indicted by the subscripts 
(XX and YY are autocorrelations). The range of E (-N, N ) 
was defined by increasing lead or lag (r) until p became zero. 
Then, p values of the correlation coefficients were calculated 
from a two-sided Student’s t distribution with the effective 
sample size (Ne) described above. We consider a > 90 % con- 
fidence level (p = 0.1) as a significant correlation, but we re- 
ported all p values when p value is less than 0. 1 allowing the 
assessments of degree of significance. Shared variance (%) 
of two time series was calculated as the square of the corre- 
lation coefficient of the two time series. 


3 Observational period (1850-2010) 

3.1 Greenland temperature anomaly 

We first investigated a reconstruction of Greenland near- 
surface air temperature in the vicinity of the Summit station 
(72°35'N, 38°27'W) during the 1850-2010 period gener- 
ated from observed temperature records and regional climate 
model output (Box et al., 2009). Uncertainties of the temper- 
ature record in earlier periods are likely larger as observa- 
tional constrains are limited (Box et al., 2009). We compared 
the record with the average temperatures of the NH (Brohan 
et al., 2006) (HadCRUT3, Fig. 1). We noted that the Sum- 
mit air temperature is highly correlated (/ =0.92, p <0.01; 
Fig. 2) with the average surface air temperature of the entire 
Greenland ice sheet (Box et al., 2009; Kobashi et al., 2011), 
indicating that the Summit temperatures reconstructed from 
ice cores are a reliable proxy for the surface temperatures of 
the entire ice sheet. Hereinafter, we call the Summit temper- 
ature time series “Greenland temperature”. Furthermore, the 
reconstructed Summit temperature for the 1855-2010 period 
is significantly correlated (r = 0.77, p <0.01) with the long 
record of observed temperatures from the Ilulissat settlement 
(Cappelen, 2011) located on the west coast of Greenland. 

The Greenland and NH average temperatures over the past 
161 yr follow similar characteristics (r = 0.66, p = 0.17; after 
linear de-trending r = 0.47, p = 0.09) but with noticeable dif- 
ferences (Box et al., 2009) as the top panel of Fig. 1 illus- 
trates. Differences in the temperature trends for Greenland 
and the NH at the decadal scale include stronger Greenland 
warming during the decade 1920-1930, decades-long cool- 
ing during the 1930s-1980s, and another warming episode 
in the 1990s and 2000s, all with similar magnitudes (Fig. 1). 
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Fig. 1. Climate records for the 1850-2010 period. Top panel: annual 
average Greenland Summit temperature (Box et al., 2009) and an- 
nual average NH temperature (HadCRUT3) (Brohan et al., 2006). 
Twenty-year root mean squared errors (RMSE) for the Summit tem- 
perature are estimated to be less than 1.6 °C (Box et al., 2009). Mid- 
dle panel: Greenland temperature anomaly (GTA[g-NH]) vs. TSI 
with background (Wang et al., 2005). Bottom panel: GTA[G-NH] 
vs. winter NAO (DJFM; CRU) (Jones et al., 1997; Osborn, 2011). 
All time series are smoothed by 3-yr RMs except the Greenland 
temperature in 51-yr RMs. Shaded areas are periods of low solar 
activity. Note that the axis for GTA is reversed. 

As a result, the decadal average temperatures in Greenland in 
the 1930s and 2000s are of similar magnitude (Fig. 1) as also 
found in changes in glacial extents in Greenland (Bjprk et al., 
2012). Unlike in Greenland, the NH temperature during the 
last decade (2001-2010) was the warmest for the observed 
period (Fig. 1). To quantitatively investigate why Greenland 
temperature deviates from the NH temperature trend, we cal- 
culated the Greenland temperature anomaly (GTA[g-nh]) by 
subtracting the standardized NH temperatures from the stan- 
dardized Greenland temperatures (Fig. 1). 



Fig. 2. The Greenland Summit surface air temperature and the 
whole Greenland ice sheet surface temperature were reconstructed 
for the last 170yr from a spatial-temporal data fusion between the 
output of a regional climate model (Polar MM5) (Box et al., 2009) 
and in situ observations (Kobashi et al., 2011). 


GTA[g-nh] = standardized Greenland temperatures 

— standardized NH temperatures 

This calculation allowed us to decompose Greenland temper- 
ature variation into three factors: background climate (NH); 
polar amplification factor (the ratio of the standard deviations 
for Greenland and NH temperatures); and regional variability 
(GTA[g— ni 1 1 )• assuming that the polar amplification factor is 
constant over the concerned period. The central Greenland 
polar amplification factor for inter-annual to multidecadal 
variation is 2.6 over the past 161 yr. 

In detail, we calculated the Greenland temperature 
anomaly (GTA[g-nh]) from the NH trend for a period of 
161 yr as follows. The variance of the inter-annual vari- 
ability of the standardized Greenland temperature is higher 
than that of the standardized NH temperature over the past 
161 yr, potentially affecting the longer-term variation of the 
GTA[g-nh]- Therefore, a 3-yr running mean (RM) filter is 
applied to the raw data for Greenland temperature and the 
NH temperature to reduce the inter-annual variability, so the 
GTA[g-nh] reflects the decadal to multidecadal temperature 
anomalies in Greenland temperature with respect to the NH 
temperature trend. After the filtering, the NH and Green- 
land temperature time series are standardized to calculate the 
GTA[g-nh]- For the standardization, we subtracted the aver- 
age value from each time series and divided the time series 
by the standard deviation of the time series; the resulting time 
series had means of zero and standard deviations of one. 

3.2 Links between NAO/ AO, GTA, and solar variability 
during the observational period 

Near-surface air temperatures in much of Greenland’s area 
exhibit a negative correlation with the NAO index especially 
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in winter (Box, 2002). As annual temperature variability 
is dominated by winter temperature variability in Green- 
land (Vinther et al., 2010), we hypothesize that the varia- 
tion in the GTA[g-nh] is linked to the NAO/AO. Indeed, the 
GTA[g-nh] is significantly negatively correlated with the av- 
erages of the winter (DJFM) NAO index (Jones et al., 1997; 
Osborn, 2011) with a correlation coefficient of r = — 0.44 
( p = 0.06) for the 1960-2007 with 3-yr RMs over the past 
48 yr as in Fig. 1. Cooling signals associated with large vol- 
canic eruptions are known to be similar to cooling signals of 
weaker solar activity, which makes it difficult to differentiate 
solar and volcanic signals (Lean and Rind, 2008; Gray et al., 
2010). However, solar signals are out-of-phase between the 
NH and Greenland temperatures (cooling in NH and warm- 
ing in Greenland, and vice versa), but volcanic signals are 
in-phase between the NH and Greenland temperatures (cool- 
ing in NH and Greenland), providing a rare opportunity to 
differentiate the two natural climate forcings. 

The inter-annual variations of the GTA[g-nh] correlates 
weakly but significantly with the 11-yr solar cycle for the 
1960-2007 period ( r = —0.29, p = 0.09), as shown in the bot- 
tom panel of Fig. 1, although the correlation is not signifi- 
cant (r = —0.27, p = 0.35) for the entire observational period 
(1851-2007). The inconsistency may relate to the empirical 
findings that climate is more sensitive to the solar forcing at 
the multidecadal-to-century time scale than shorter variations 
of solar activity such as the 11-yr solar cycle (Waple et al., 
2002; Gray et al., 2010). In the multidecadal-to-centennial 
time scale, the GTA[g-nh] indicates an increasing trend un- 
til the 1930s. Note the GTA[g-nh] vertical axis is reversed 
in Fig. 1. The GTA[g-nh] thereafter showed a decreasing 
trend, indicating that Greenland temperature decreased rel- 
ative to that of the NH until ~ 1990. Then, the GTA[g-nh] 
increased to the present associated with more rapid warming 
in Greenland than that in the NH (Fig. 1). The multidecadal 
trend of the GTA[g-nh] is similar to the long-term variation 
of solar activity over the past 160 yr (Fig. 1). However, the 
observational records are too short to support a causal rela- 
tionship between the solar variability and the GTA[g-nh] in 
the multidecadal-to-centennial time scale. 


4 The past 800 yr of Greenland temperature anomaly 
(GTA[g-nh]) 

To evaluate the robustness of the apparent solar influence on 
the multidecadal-to-centennial GTA[g-nh] variation, we ex- 
amined records from the past 800 yr (Fig. 3). For the NH 
temperature, two proxy-based reconstructions by Moberg et 
al. (2005) and Mann et al. (2009) were used as in Fig. 3. For 
Greenland temperature, we used a recently developed record 
of the temperature at Summit that was reconstructed from 
air trapped in the GISP2 ice core (Kobashi et al., 2010). The 
new temperature record addresses deficiencies in the existing 
Greenland temperature proxies, such as <5 1 s Oj ce and borehole 


thermometry (Kobashi et al., 2011), providing decadal av- 
erages owing to gas diffusion in the firn (unconsolidated 
snow) layer. Thus, this new temperature record is not af- 
fected by large seasonal temperature swings (Kobashi et al., 
2010). Owing to gas and heat diffusion in the firn layer 
(porous unconsolidated snow), the temperature record rep- 
resents decadal averages (Kobashi et al., 2010). For the as- 
sessments of solar influence on <$ 18 Oi ce , we refer to Kobashi 
et al. (2012). 

The NH and Greenland temperatures exhibit similar 
trends, with the characteristics of the Medieval Warm Period 
in 1200-1450, the Little Ice Age in 1450-1850, and recent 
warming (Fig. 3). We then derived two 800-yr GTAs[g-nh] 
by the same method used for the observational records, but 
fully taking into account uncertainties in temperature recon- 
structions (Fig. 3). For the data over the past 800 yr, we ini- 
tially applied 21 RMs to the NH and Greenland temperatures 
prior to standardization to illustrate multidecadal to centen- 
nial variability in the GTA[g-nh]- We standardized the time 
series by subtracting the means and dividing the time series 
by the standard deviations, so the resulting time series have 
means of zero and standard deviations of 1 . Ninety-five per- 
cent confidence intervals for the NH (Mann et al., 2009) and 
Greenland (Kobashi et al., 2010) temperatures with 21 RMs 
were estimated by Monte Carlo simulation as follows. Gaus- 
sian white noise signals with a standard deviation of 1 and 
a mean of zero are generated and then multiplied by the re- 
ported la values for the NH (Mann et al., 2009) and Green- 
land (Kobashi et al., 2010) temperature data for each year. 
These modified signals were then added to the original NH 
and Greenland temperature time series. By repeating this cal- 
culation 1000 times, 1000 synthetic NH and Greenland tem- 
perature time series were generated. For each time series, 
21 RMs were applied. The ler values of the NH and Green- 
land temperature time series with the 21 RMs were then es- 
timated as standard deviations from the 1000 synthetic time 
series. 

To address the uncertainty in the NH temperatures in 
the 21 RMs of Moberg et al. (2005), we directly used 
the uncertainty (uncertainty A) due to the variance among 
low-resolution proxies, but we did not include the uncer- 
tainty (uncertainty B) in the scaling factor or the constant 
adjustment term (uncertainty C), because these uncertain- 
ties were not relevant for standardizing time series. The 
estimated la values of the NH and Greenland tempera- 
tures with 21 RMs were scaled for the standardized NH 
and Greenland temperatures in 21 RMs using the same 
values for standardizing the NH and Greenland tempera- 
tures with 21 RMs. The uncertainties in the GTA[g-nh] 
were derived from the equation for the propagation of 
error: aGTA = y 7 (a stn. NH) 2 + (a stn. Greenland) 2 , where 
<tGTA[g-N[I|. a stn. NH, and a stn. Greenland are the stan- 
dard deviations of the GTAjg-nh]- the standardized NH 
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Fig. 3. Climate records for the 1200-2010 period, (a) Greenland temperature (blue) (Kobashi et al., 2010) and the NH temperature from 
Moberg et al. (2005) (red) and Mann et al. (2009) (green) in 21 RMs with 95 % confidence intervals (2cr). (b) The blue line is the standardized 
Greenland temperature (Kobashi et al., 2010). The red and green lines are the standardized NH temperatures (Moberg et al., 2005; Mann 
et al., 2009). (c) Greenland temperature anomalies (GTA[q_nh]) calculated from Greenland temperature and the NH temperature, as shown 
in panel (b). The blue and red lines are smoothed by 21 RMs with 95 % confidence intervals (2a) and by 51 RMs, respectively. The area 
between the red line and zero is colored blue for the region above zero and red for the region below zero. The black line is the GTA[q_nh] 
calculated from the observed temperatures (5 1 RMs ) in Fig. 1 . Transparent light blue shading indicates areas below the observed GTA[q_nh] 
(black line) and above zero. Values of — 1 .27 and —0.49 are added to the observed GTAs[g_nh] with the NHs by Moberg et al. and Mann 
et al., respectively, to account for the differences in the average GTAs[q_nh] between the proxy and the observations for the overlapping 
period. Note that the y-axis of the GTA[g_nh] is reversed to allow for comparison with the solar variability, (d) The purple line represents 
the solar forcing from Bard et al. (2000) and Crowley (2000) spliced into that from Lean et al. (1995), and the green line is from Delaygue 
and Bard (2009) spliced into that from Wang et al. (2005) The areas are colored red when the values are above zero (stronger sun) and blue 
when the values are below zero (weaker sun). The standardized time series of the solar forcings is shown. The grey areas are the periods of 
weaker suns used for the statistical tests, (e) The volcanic forcing for the NH from Gao et al. (2008) is plotted as a reference. 
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temperature, and the standardized Greenland temperature, 
respectively. 

The two GTAs[g-nh] exhibit similar characteristics over 
the past 800 yr (r = 0.79, p = 0.02), partly because two NH 
temperature records closely resemble each other, indicat- 
ing that the temperature signal is a robust feature of the 
GTA[g-nh]- The GTAs[g-nh] explains 31-35 % of the vari- 
ation of Greenland temperature over the past 800 yr in the 
multidecadal to centennial time scale. The central Green- 
land polar amplification factor (=Greenland/NH) over the 
past 800 yr ranges from 3.3 for the NH temperature from 
Moberg et al. (2005) to 4.2 for the NH temperature from 
Mann et al. (2009). The GTAs[g-nh] show more negative 
values (with Greenland being cooler than the NH) before 
1400 and after 1700 and positive values (with Greenland be- 
ing warmer than the NH) between 1400 and after 1700, indi- 
cating that there was a more positive NAO-like pattern before 
1400 and after 1700, and a more negative pattern in the inter- 
vening period. The persistent positive NAO during the Me- 
dieval Warm Period is consistent with the results of a recent 
study (Trouet et al., 2009). 


5 Solar influence on Greenland temperature anomalies 
(GTA[g-nh]) 

To investigate the relationship between the GTAjg-nh] and 
changes in solar activity over the past 800 yr, we used two 
solar forcing reconstructions based on 10 Be in ice cores from 
Antarctica (Bard et al., 1997; Crowley, 2000). Cross wavelet 
analyses (Grinsted et al., 2004) in Fig. 4 revealed that Green- 
land temperature has a weak anti-phase correlation with the 
total solar irradiance (TSI) around the period of ~ 128 yr. In 
contrast, the NH temperature generally exhibits a weak in- 
phase relationship with the TSI over the past 800 yr (Fig. 4). 
As a result, solar signals are enhanced in the GTA[g-nh] and 
exhibit significant anti-phase with the TSI, especially around 
the period of ~ 128 yr (Fig. 4). 

To test the null hypothesis that solar variation did not cause 
changes in GTAs[g-nh] over the past 800 yr, we conducted a 
statistical test “re-sampling” (Ito and Minobe, 2010). We de- 
termined the periods of stronger and weaker solar output us- 
ing two times series of solar forcing, as shown in Fig. 3. The 
two time series of solar forcing were smoothed by 21 RMs, 
and then standardized. The two resulting time series were av- 
eraged to generate a representative time series of solar forc- 
ing. We simply defined the periods of stronger solar output as 
periods when the standardized solar forcing was positive, and 
the periods of weaker solar output were defined as the periods 
when the standardized solar forcing was negative (Table 1). 
In total, there were 384 yr of sttonger solar output, and 417 yr 
of weaker solar output. 

Then, the averages of the GTA[g-nh] corresponding to the 
periods of “stronger” and “weaker” solar output were cal- 
culated. The averages of the GTA[g-nh] for stronger and 


Table 1 . The periods of stronger and weaker solar output over the 
past 800 yr. 


# 

Solar 

output 

Top 

Bottom 

Duration 

1 

Strong 

2000 

1826 

175 

2 

Weak 

1825 

1800 

26 

3 

Strong 

1799 

1740 

60 

4 

Weak 

1739 

1625 

115 

5 

Strong 

1624 

1597 

28 

6 

Weak 

1596 

1396 

201 

7 

Strong 

1395 

1345 

51 

8 

Weak 

1344 

1270 

75 

9 

Strong 

1269 

1200 

70 


weaker solar output were subtracted to obtain the average 
difference in the GTA[g-nh] between the periods of stronger 
and weaker output. Second, the order of the 1 1 durations in 
Table 1 was randomly shuffled. Based on the new order of the 
durations, 1 1 periods were calculated from 2000 to 1200CE. 
Periods of stronger and weaker solar output were defined as 
the odd and even numbers, respectively, in the new order. In 
addition, the assignments of odd and even numbers to pe- 
riods of stronger and weaker solar output were randomly re- 
placed. A synthetic GTA[g-nh] is generated by adding Gaus- 
sian white noise (with a zero mean and a standard deviation 
of 1) multiplied by the uncertainties (ler) of the previously 
calculated GTAs[g-nh]- The averages of the GTAs[g-nh] 
from the periods of stronger and weaker solar output were 
then calculated, and the difference between these two aver- 
ages was calculated. This process was repeated 10000 times 
(1000 times for maps) to generate a set of 10000 samples 
(1000 samples for maps) for the differences. Percentile ranks 
were calculated and used to test the null hypothesis with one- 
sided tests as we are only interested in Greenland cooling 
(warming) not warming (cooling) during stronger (weaker) 
sun, respectively, to determine the significance levels of the 
observed differences in the GTA[g-nh] between periods of 
stronger and weaker solar output. We found that, for the 
GTA[g-nh] with respect to the NH temperature of Moberg 
et al. (2005), the null hypothesis was rejected at the 99.9 % 
significance level; for the GTAjg-nh] with the NH tempera- 
ture of Mann et al. (2009), the null hypothesis was rejected 
at the 98 % confidence level. Therefore, we conclude that so- 
lar variability has caused changes in the GTA[g-nh] over the 
past 800 yr. 

One caveat for the calculation of the GTA[g_nh] from 
Greenland temperature and the NH temperature is that the 
correlation of the GTA[g-nh] with solar activity may arise 
entirely from the correlation between the NH temperatures 
and solar activity. To test the null hypothesis that the correla- 
tions between the GTA[g-nh] and solar variability originate 
completely from the positive correlation between the NH 
temperature and solar variation, we conducted a statistical 
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Fig. 4. Cross-wavelet transform (XWT) (Grinsted et al., 2004) of Greenland temperature, the NH temperature, the GTA[g-nH]> and the 
NAO, all vs. TSI, for 1200-2010. (a) Greenland temperature (Kobashi et al., 2010) vs. TSI (Lean et al., 1995; Bard et al., 2000; Crowley, 
2000) for 1200-2010; (b) the NH temperature (Moberg et al., 2005) vs. TSI (Lean et al., 1995; Bard et al., 2000; Crowley, 2000) for 1200- 
2010; and (c) the GTA[G-NH] with the NH temperature from Moberg et al. (2005) vs. TSI (Lean et al., 1995; Bard et al., 2000; Crowley, 
2000) for 1200-2010; (d) the NH temperature (Mann et al., 2009) vs. TSI (Lean et al., 1995; Bard et al., 2000; Crowley, 2000); (e) the 
GTA[q_nh] (Moberg et al., 2005) vs. TSI (Lean et al., 1995; Bard et al., 2000; Crowley, 2000). Annual resolution data are used for the 
analyses except the GTAsjg-NH]- The area between the grey lines represents the 95 % confidence level relative to the noise (red). Arrows 
indicate relative phase relationships, with in-phase relationships denoted by right-facing arrows; arrows pointing straight down indicate that 
the phase of the TSI leads that of Greenland temperature, the NH temperature, or the GTA[G-NH] by 90°. Shaded areas represent the cone 
of influence where edge effects may distort the results. 


test by decomposing the standardized Greenland temperature We note that 7b = GTA[g-nh] ■ Therefore, we can write it as 
(f G ) and the standardized NH temperature (7 nh) as follows: the following: 


Ta = (j'G + Tnh) j \/2 

7b = (lb - 7 nh) j V2 . 


T G = (T a + 7b) /V 2 


Tnh = (7a - T b )/V2. 
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Table 2. Degrees (9) of the maximum negative correlation coefficients, r$(9), in various settings. The 95 % confidence intervals (2a) are 
shown for the past 800 yr with the T^h from Moberg et al. (2005) and Mann et al. (2009), including uncertainties in Greenland temperature 
and the NH temperature. For the past 800 yr, two sets of solar forcing were used: (a) the dataset from Bard (Bard et al., 2000) and Crowley 
(Crowley, 2000) spliced into that of Lean et al. (Lean et al., 1995), and (b) the dataset from Delaygue and Bard (Delaygue and Bard, 2009) 
spliced into that of Wang et al. (Wang et al., 2005), as shown in Fig. 3. 


Datasets 


Degrees (9) of maximum negative correlation, r s (9) 


Past 161 yr (observational) 

Past 1000 yr; by Moberg et al. (2005) 
Past 1000 yr; by Mann et al. (2009) 
Simulated past lOOOyr; GISS (ensemble) 


129° 

(a) 11 1.7 ± 1.0°; (b) 111.7 ± 1.1° 
(a) 1 13.8 ± 0.8°; (b) 113.9 ±0.9° 
114° 



B.C.E. 


Years 


C.E. 


Fig. 5. Greenland temperature and solar activity over the past 4000 yr. Both time series are linearly de-trended and smoothed with 21-yr 
RMs. Greenland temperature (Kobashi et al., 201 1) and solar activity (Steinhilber et al., 2009) have weak but significant negative correlation 
(r = —0.18; p = 0.09). 21-yr RMs were chosen as other smoothing with various periods showed lower confidence levels although they had in 
some cases higher correlation coefficients (see panel in left bottom). Note that the axis for the TSI is reversed. 


Here, we define T c (0) as 
T c (6 ) = T. d cos 6 + 7t, sin 9. 
Clearly, 


7’;, 


= 0 ° 


T c {0) = 


f G ; 0 = 45° 

7b = GTA; 0 = 90° 

— 7nh; 0 = 135° 

-r a ; 0 = 180° 


correlation occurs around 110-114° (Table 2). The results 
of the GISS simulations yielded a value similar to the ob- 
served value (Table 2). Therefore, we concluded that when 
solar activity is stronger (weaker), GTA[g-nh] becomes 
lower (higher) due to higher (lower) NH temperatures and 
lower (higher) Greenland temperatures, respectively with the 
99.0% significance level. The finding is also consistent with 
the fact that the Greenland temperature has a weak but sig- 
nificant negative correlation with solar activity over the past 
4000 yr as shown in Fig. 5 and Kobashi et al. (2012). 


We then examined the correlation (r s (0)) between changes 
in solar forcing and T c (0 ) as r s (0) = (solar forcing, T c (0)) 
under variable 0. If the negative correlation between the 
GTA[g-nh] and solar forcings originates completely from 
the positive correlation between 7 \jh and solar forcing, the 
maximum negative correlation of T c (0) should occur at 
0 > 135°. In contrast, if Tq contains anti -phase components 
with solar forcing and thus contributes to the negative corre- 
lation between the GTA[g-nh] and solar forcing, the maxi- 
mum negative correlation of r s (0) should occur in the range 
45° <6 < 135°. For the past 800 yr, the maximum negative 


6 Relationship between NH and Greenland 
temperatures with solar variability 

To further investigate the relationship between NH and 
Greenland temperature, we calculated another principal 
component of the two time series, “GTP[g-nh]” (Green- 
land temperature plus = standardized Greenland tempera- 
ture + standardized NH temperature). In contrast to the GTA 
that illustrates out-of-phase climate signal (solar) of the 
two time series, GTP[g-nh] can be used to investigate in- 
phase climate signals in the two time series. Therefore, the 
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Years (C.E.) 


Fig. 6. Results of a simple numerical modeling. Top panel: 
GTP[ G+N h] and millennial trend from polynomial regression of 
GTP [G +NH] with the NH temperature by Mann et al. (2009). The 
volcanic forcing for the NH from Gao et al. (2008) is plotted as a 
reference. Shaded areas are low GTP periods in the multidecadal- 
to-centennial time scale. Note that many shaded areas coincide 
with large volcanic eruptions. Middle panel: synthetic and observed 
Greenland temperature (Kobashi et al., 2010) and NH temperature 
(Mann et al., 2009). Bottom panel: synthetic GTAjg-nhr observed 
GTAs[ G _nh] with the NH temperatures by Mann et al. (2009) and 
Moberg et al. (2005), and standardized solar forcing (solar 2 in 
Fig. 3). 


GTP[g-nh] variation (Fig. 6, top panel) in the multidecadal- 
to-centennial time scale likely reflects changes in volcanic 
forcing, increasing greenhouse gas forcing, and possibly in- 
ternal variation of climate. Millennial-scale variation (Fig. 6, 
top panel) of the GTP[g-nh] peaking in the 1300s and 1700s 
may reflect non-linear climate responses (e.g., growth of 
glaciers and ice sheets) to solar and/or volcanic forcing. 

To illustrate the relationship between the GTA[g-nii|. 
GTP|g+mi|, and solar activity, we conducted a simple nu- 
merical experiment. To start with a synthetic hemispheric 


temperature without solar signals, a trend equation was de- 
rived from polynomial regression of the GTP[g+nh] by the 
NH temperatures by Mann et al. (2009) and the Greenland 
temperature (Kobashi et al., 2010) (Fig. 6). These values of 
trend for a period of 1200-2000 were multiplied by 1.3. 
Then, the random numbers that were greater than or equal 
to 0 and less than 1 were generated, and multiplied by 8. 
Then, these random numbers were added to the trend to pro- 
duce seed time series. Synthetic NH temperature was pro- 
duced by adding the seed to the standardized solar signal 
(solar 2 in Fig. 3) multiplied by two. Then, the synthetic NH 
temperature was smoothed by 11-yrRMs, and standardized. 
Synthetic Greenland temperature was generated by adding 
negative solar signals (multiplied by 3) to the synthetic NH 
temperatures. We note that all the coefficients were obtained 
to match the observed temperature records. 

Then, the synthetic Greenland temperature was smoothed 
by 11 -yr RMs, and standardized. Synthetic GTA[g-nh] was 
calculated by subtracting the standardized synthetic NH tem- 
perature from the standardized synthetic Greenland temper- 
ature. In addition, synthetic GTP[g+nh] was calculated by 
adding the standardized synthetic NH temperature and the 
standardized synthetic Greenland temperature. The fact that 
residuals between the model results and the observations (20- 
values are 1.6 for the standardized Greenland temperatures 
and 1.2 for the standardized NH temperatures as in Fig. 6, 
middle panel) are smaller than the long-term changes (3 for 
the standardized Greenland temperatures; 2.7 for the stan- 
dardized NH temperatures as in Fig. 6, middle panel) sup- 
ports the validity of the model. 

We found that the synthetic GTA[g-nh] captures the orig- 
inal solar signal (r = — 0.83) with slightly different features. 
The synthetic GTA[g-nh] shows slightly higher correlation 
(r = 0.59) with the observed GTA[g-nhj than the correlation 
(r = —0.55) with the original solar signal (Fig. 6), indicat- 
ing that the GTA[ G -nh] calculation introduces a minor alter- 
ation on solar signal as in Fig. 6. The synthetic GTP[ G +nh] 
can fully recover the original millennial trend (r > 0.99). In 
addition, the synthetic GTA[g-nh] calculated with synthetic 
NH temperature without solar signal (the trend and noises 
only) and synthetic Greenland temperature (with negative so- 
lar signal) is more strongly correlated with the original solar 
signals, indicating that the observed correlation between the 
solar forcing and the GTA[g-nh] is not an artifact of solar 
signal in the NH temperature. 


7 Climate model experiments 

We examined the climate model responses of the 
GTAs[ G -nh] in a coupled atmosphere-ocean climate 
model, the NASA Goddard Institute for Space Studies-ER 
(GISS-ER) simulation (Schmidt et al., 2006; Mann et al., 
2009). The GISS model used in this study is the mode used 
in the CMIP3 simulations performed for the IPCC AR4 
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(Schmidt et al., 2006). The model was run at 4 x 5 degrees 
of horizontal resolution with 23 vertical layers (many of the 
CMIP3 runs used a 20-layer version without gravity-wave 
drag but were otherwise identical). We performed a long 
control run to establish stable initial conditions, from which 
six transient runs were branched off, each extending from 
850 to 1900 CE. Solar forcing was applied across the 
ultraviolet, visible and infrared spectra based on scaling by 
wavelength versus total irradiance using the relationships as 
observed in modern satellite data. The total irradiance over 
time was based on the time series of 10 Be data from a South 
Pole ice core (Bard et al., 2007) and by taking into account 
a small long-term geomagnetic modulation (Korte and Con- 
stable, 2005) and a polar enhancement factor. The amplitude 
was scaled to provide a top-of-the-atmosphere irradiance 
change of 1.1 W m -2 between the Maunder Minimum and 
the late 20th century, corresponding to a tropopause radiative 
forcing of 0.38 W/m 2 as in Shindell et al. (2006) and Wang 
et al. (2005). The model also includes the ozone photo- 
chemical response to solar irradiance variations, which is 
parameterized from the results of prior GISS modeling using 
a full atmospheric chemistry simulation (Shindell et al., 
2006) for computational efficiency. The results from these 
millennial transient simulations and those from a second set, 
including doubled solar forcing magnitude performed to test 
the sensitivity of the climate response, were also presented 
by Mann et al. (2009) and by Pechony and Shindell (2010). 

The GTAs[g-nh] in the model were calculated using 
the average temperature of southern Greenland including 
the GISP2 site as the northern limit (72° 36' N, 38°30' W; 
3200 m a.s.l.). The ensemble averages of the six GTAs[g-nh] 
calculated with GISS-ER model were compared with the 
changes in solar irradiance (Fig. 7), and a similar corre- 
lation (r = — 0.47, ranging from r = 0.11 to —0.75 in in- 
dividual runs) to the observation was found between the 
GTA[g-nh] and solar variability in the multidecadal-to- 
centennial timescale (e.g., 51-yr RMs; Fig. 7). The re- 
sampling test indicated a significant solar influence on the 
GTA[g-nh] variation in the GISS-ER model in a 97 % confi- 
dence level. We also tested the null hypothesis that the corre- 
lation between the GTA[g-nh] and solar forcing in the model 
is a result of correlations between the NH temperature and 
the TSI. The null hypothesis was rejected for the GTA[g-nh] 
in the GISS-ER simulation, indicating Greenland tempera- 
ture negatively responded to changes in solar activity. There- 
fore, the observed GTA[g-nhj variations, including the neg- 
ative correlation between Greenland temperature and solar 
forcing, are physically explainable. This result is consistent 
with the results of earlier studies showing that the GISS sim- 
ulation with integrated ozone feedback processes appears to 
properly reproduce many aspects of solar-induced NAO/ AO- 
like dynamics (Mann et al., 2009). 



Fig. 7. Standardized solar variations (top panel) and GTAsjg-nH] 
in the GISS-ER simulations (bottom panel). The solar variations 
and shaded areas are the same as in Fig. 3. Thin lines in the 
GTAs[G-NH] represent each individual run, and the thick black line 
is the ensemble averages. Note that the y-axes for the GTAs[G-NH] 
are reversed. 

8 Dynamic causes of solar-induced GTA[q-nh] 
variation 

To investigate the spatial extent and mechanisms of solar- 
induced GTAjg-nh] variations, we plotted the composite- 
mean temperature differences of periods of stronger and 
weaker solar output, as defined in Fig. 3, for the model out- 
puts and the grid proxy reconstruction of global temperatures 
(Mann et al., 2009) in Fig. 8. Importantly, the spatial features 
in the proxy map resemble the surface air temperature re- 
sponses to the 11-yr solar cycles for the 1889-2006 period 
(Lean and Rind, 2008). The locations and extents of small 
areas of cooling in the northern North Atlantic are consis- 
tent in the two maps (Fig. 8). The cooling signal cannot be a 
result of the contamination of volcanic signals as large vol- 
canic events were observed to cause positive NAOs and thus 
an in-phase relationship between the NH and Greenland tem- 
peratures. In addition, note that the temperatures in southern 
Greenland in the GISS-ER show cooler areas. 

Earlier studies using the GISS-ER simulation and obser- 
vations have shown that the NAO/ AO-like temperature pat- 
tern is related to changes in large-scale dynamic atmospheric 
circulation (Lean and Rind, 2008) and polar night jets re- 
sulting primarily from solar spectral (UV) changes (Kodera 
and Kuroda, 2002; Kuroda and Kodera, 2002). Composite 
differences in the sea level pressure (SLP) in the GISS-ER 
model indicate significant pressure increases over both the 
Atlantic and Pacific sectors at mid-latitudes and decreases 
over the majority of the northern high latitudes (Fig. 9). Al- 
though the SLP anomalies show a more complex structure 
than the NAO and AO patterns and the signals are rather 
weak, these changes do project onto both of those variability 
modes and lead to altered flow from North America towards 


Clim. Past, 9, 583-596, 2013 


www.clim-past.net/9/583/2013/ 


T. Kobashi et al.: On the origin of multidecadal to centennial Greenland temperature anomalies 


593 



-0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2 


Fig. 8. Composite annual mean temperature differences (stronger 
solar output - weaker solar output) as defined in Fig. 3 from the 
model output (GISS-ER) and the proxy data (Mann et al., 2009). 
Only grids with 95 % significance levels, calculated by re-sampling 
with two-sided tests as we do not have a priori knowledge on 
warmer or cooler responses to solar variations except in Greenland, 
are shown in color. For the GISS-ER (ensemble mean result), the 
solar forcing from the Maunder Minimum to the late 20th century 
is 0.38 Wm“ 2 . 

Greenland and the North Atlantic, similar to that induced by 
NAO/AO anomalies (Fig. 9). As a result, Greenland experi- 
ences stronger northerly winds and resultant cooler tempera- 
tures. Earlier studies with GISS-ER and a previous version 
of GISS have shown that the solar-induced NAO/AO-like 
patterns are amplified by the interactive ozone processes in 
the atmosphere, which is called “top-down process” (Gray 
et al., 2010) as described previously (Shindell et al., 2001; 
Mann et al., 2009), and may also be linked to surface changes 
for instance in tropical and subtropical moisture convection 
processes through surface heating by variable solar activity 
(Shindell et al., 2001; Kodera, 2002; Kodera and Kuroda, 
2002 ). 

Regarding the observed centennial to multidecadal varia- 
tions in GTAs[g-nhf the mechanisms may involve oceanic 
processes (Cubasch et al., 1997; Waple et al., 2002) such 
as the Atlantic meridional overturning circulation (AMOC). 
To further investigate the dynamic causes of the GTA[g-nh] 
variations, we investigated changes in the AMOC. The 
GISS-ER simulation exhibits a significant reduction of 0.3 



- 0.2 - 0.15 - 0.1 - 0.05 0 0.05 0.1 0.15 0.2 

Fig. 9. Composite SLP differences in the GISS simulation (using 
the same data as in Fig. 8 except for the SLP). Grid data with 90 % 
confidence levels, calculated by re-sampling (two-sided test), are 
shown in color. 


±0.1 Sv from the mean value of 21.1 Sv in the AMOC 
(evaluated across 48° N at 900 m depth) as a difference be- 
tween periods of stronger and weaker solar activity as defined 
above, thus indicating slowed heat transport from the South 
to the North Atlantic, which may explain part of the anti- 
phase relationship between the NH and the southern Green- 
land temperatures. The same AMOC processes are often ob- 
served in model projections of future warming (Meehl et al., 
2007). Owing to reduced AMOC in many models, northern 
North Atlantic areas including northern Europe are projected 
to experience less warming. 

To further investigate the relationship between AMOC and 
GTA, we looked into a reconstructed North Atlantic average 
temperature record (Mann et al., 2009) over the past 800 yr 
(Fig. 10). Greenland Summit temperatures exhibit significant 
correlation with average North Atlantic temperature over the 
observational period as in Fig. 4 in Kobashi et al. (2012) 
with a correlation coefficient of /• = 0.58 (p = 0.03) in the 
annual time scale, and over the past 800 yr (Fig. 10) with 
r = 0.60 ( p = 0.09) in the multidecadal time scale. Multi- 
decadal variation of average North Atlantic temperature, of- 
ten called Atlantic multidecadal oscillation (AMO), has been 
linked with changes in AMOC variation and resultant inter- 
hemispheric heat transport in North Atlantic (Delworth and 
Mann, 2000; Knight et al., 2005; Ottera et al., 2010). We cal- 
culated North Atlantic temperature anomalies (NATAs) sim- 
ilar to GTA. It is found that NATAs possess similar charac- 
teristics with GTA (r = 0.5, p = 0.1 8, Fig. 10, bottom panel), 
indicating North Atlantic temperature was relatively warmer 
during the 15th century (at the time of lower solar activity 
during the Sporer Minimum) and cooler during the past few 
centuries (higher solar activity) than the NH temperature. We 
note that volcanic events may also have played an important 
role for the variation of NATA (Ottera et al., 2010), as more 
volcanic eruptions occurred during the lower solar activity 
over the past millennium. We speculate that, during lower 
(higher) solar activity, AMOC was stronger (weaker), and so 
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Fig. 10. Greenland and North Atlantic temperature anomalies over 
the past 800 yr. Top panel: Greenland temperature (Kobashi et al., 
2010) and North Atlantic average temperature anomaly (Mann et 
al., 2009) in 21-yr RMs. Middle panel: standardized North Atlantic 
and NH average temperatures (Mann et al., 2009) in 21-yr RMs. 
Bottom panel: GTA (Mann NH) in 21-yr RMs is from Fig. 3c. North 
Atlantic temperature anomaly (NATA) was calculated by subtract- 
ing standardized NH temperature from standardized North Atlantic 
average temperature as in the middle panel. It is noted that Green- 
land temperature is significantly correlated with North Atlantic av- 
erage temperature ( r =0.60, p = 0.09 in 21-yr RMs), while it is not 
with NH average temperature (r = 0.43, p = 0.21 in 21-yr RMs). 
Ninety-five percent confidence intervals (2a) are estimated as for 
Fig. 3. Shaded areas represent periods of low solar activity same as 
in Fig. 3. 

temperature was relatively higher (lower) in North Atlantic 
than the NH temperature over the past 800 yr. Therefore, the 
observed negative responses to solar variation in the northern 
North Atlantic (Greenland) suggest that oceanic responses 
to changes in solar outputs played a role in the GTA[g-nh] 
variations. 


9 Conclusions 

Over the past 800 yr, Greenland temperature variability was 
dominated by background climate changes at the northern 
hemispheric scale, and these changes are known to be largely 
the result of solar and volcanic forcing prior to 1850 (Crow- 
ley, 2000). However, we find that the deviations in Greenland 
temperatures from that of the NH are also strongly influ- 
enced by atmospheric dynamics such as those produced by 
the NAO/AO, in which longer-term variations are induced 
by atmospheric circulation responses to variations in solar 
outputs, with an additional contribution from the oceanic re- 
sponse. This finding indicates that the Greenland temperature 
minimum over the past 800 yr occurred in the early 18th cen- 
tury, in contrast to the minimum for the NH temperature in 
the 15— 16th century (Figs. 3 and 6) - likely because of in- 
creased solar forcing during the 18th century. In addition, 
during the latter half of the 20th century, at time of high solar 
activity, Greenland temperature exhibited less pronounced 
warming than the NH. If future solar activity were to be low 
compared with that of the past lOOyr (Livingston and Penn, 
2009), it may reveal a hidden warming of Greenland caused 
by increasing greenhouse gases during the latter half of the 
20th century. 
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